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Abstract 

Differential cross sections of neutral pions emitted in ^'^^Ta+^^^Au collisions at a beam energy 
of 39. 5A MeV have been measured with the photon spectrometer TAPS. The kinetic energy and 
transverse momentum spectra of neutral pions cannot be properly described in the framework of 
the thermal model, nor when the reabsorption of pions is accounted for in a phenomenological 
model. However, high energy and high momentum tails of the pion spectra can be well fitted 
through thermal distributions with unexpectedly soft temperature parameters below 10 MeV. 
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I. INTRODUCTION 



Neutral pion emission is a very sensitive probe of the processes taking place in the inter- 
action zone of two nuclei colliding at beam energies below lOOA MeV. This energy is well 
below the free nucleon-nucleon (NN) vr*^ threshold of 280 MeV. In consequence, the energy of 
produced mesons is limited, even with the help of cooperative effect of Fermi motion. This 
energy limitation allows for the experimental study of the energy spectrum and angular 
distribution of 7r° mesons to a large extent. In addition, the dominant two-photon decay 
channel of this short lived meson makes possible the detection of neutral pions also at rest. 
The number of detected neutral pions, superior to the previous measurement [l, 2], allows 
us to make more definitive conclusions about their spectral properties. 

With decreasing beam energy per nucleon, the particle production process increasingly 
relies on the internal motion of nucleons and their correlations (cooperative effects) as well 
as on the dynamical evolution of the colliding system. The particles produced in these 
unfavourable conditions carry important information about the collision process, as they 
remove a significant fraction of the available energy. Various models have attempted to 
describe these processes. Fermi motion, which plays a crucial role in subthreshold particle 
production, is, of course, incorporated in modern transport model calculations {e.g. j^). 
These models include baryonic resonances, which serve as a temporary energy storage en- 
abling the production of mesons. However, while the gross properties of subthreshold particle 
production are relatively well explained by these models, the description of spectral proper- 
ties of neutral pions does not always provide satisfactory results. While for Ar+Al collisions 
at 95A MeV a rather consistent description of 7r° emission was found [4] , the yield of neutral 
pions from Kr-|-Ni collisions at 60A MeV was strongly underpredicted in the high energy 
part of the spectrum Q] . Phenomenological models of first-chance nucleon-nucleon collisions 
have also been used to describe the energy scaling of the cross sections {e.g. ol). In another 
extreme, purely thermal models of particle emission have been applied to describe the inter- 
action zone of two colliding nuclei [7|. Strong reabsorption processes of pions in the nuclear 
medium have been incorporated in both the transport and thermal models fl^, ^• 

The production of neutral pions was studied with the TAPS spectrometer 91] for a nearly 
symmetric heavy system -'^'^-'^Ta+^^''Au at a beam energy of 39. 5A MeV, i.e. at 1/7 of the free 
NN 7r° production threshold energy. Particular attention was paid to removing all effects that 
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could obscure the experimental data for such rarely produced particles. The contribution 
from cosmic radiation was reduced to below 1% of the final neutral pion sample. 



II. EXPERIMENTAL SETUP 

In the experiment performed at GANIL ^^^Ta^'^+ ions were accelerated to 39. 5 A MeV 
beam energy and delivered in bunches with a rate of 9.055 MHz. The nominal beam intensity 
was 10 enA, ca. 160 ions per pulse. The Ta ions impinged on a 19.3 mg/cm^ thick ^^^Au 
target. The interaction probability per average beam pulse was about 7%. 

Photon pairs originating from tt" decays were detected by the TAPS electromagnetic 
calorimeter jo], arranged in 6 blocks, each composed of 64 hexagonal modules in an 8x8 
honeycomb structure. The configuration of the experimental setup was shown in Fig. 1 of 
Ref. |10|] (in our measurement the inner detectors Dwarf Ball and SSD were removed). The 
blocks covered a wide range of polar angles (between 45° and 170°), making up ca. 20% of 
the full solid angle. Each module consisted of a 25 cm-long BaF2 crystal (corresponding 
to 12 radiation lengths) of hexagonal cross section with an inscribed radius of 2.95 cm. 
A charged particle veto detector (CPV) made of plastic scintillator NE102A was mounted 
in front of each module. Events with neutral pion candidates were selected by a trigger, 
requiring activation of at least two modules by neutral particles, in different TAPS blocks. 
An activation was defined by deposing energy of about 12 MeV, a value not affecting the 
measurement of energy distribution of photons originating from neutral pion decays. 

III. DATA ANALYSIS 

A single TAPS module measures the time, energy and electric charge state (neutral or 
charged) of the incoming particle. The scintillation light response of BaF2 has two com- 
ponents, whose relative intensity differs for electrons (and photons) compared to heavily- 



ionizing particles like protons and deuterons ll|]. This feature, combined with the Time of 



Flight (ToF) infornaation and the charge veto signal, delivers nearly unambiguous identifica- 
tion of a photon [12]. This particle, traversing the BaF2 module, induces an electromagnetic 
cascade which usually spreads into adjacent modules, creating a cluster of active modules. 



The properties of a particle were reconstructed from a cluster using the ROSEBUD 
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FIG. 1: Top: invariant mass distribution of photon pairs shown with fit of an asymmetric Gaussian 
function to the vr*^ peak and an exponential background. Hatched histogram: the vr*^ candidates 
after the kinematical fit selection; dashed histogram: contamination from cosmic radiation. Bot- 
tom: Phase space distribution of the reconstructed vr^ in the px vs yiab plane; dashed curve for the 
average momentum of the vr^ obtained from a Boltzmann distribution characterized by T=15 MeV; 
the arrow indicates the Vcf^ (see text). 

analysis package. If a particle hits the border of a TAPS block, part of the induced cascade 
can escape the detector. Such a cluster was rejected by requiring that the maximum energy 



deposition can not be in a border module (c.f. [ISj]). Cosmic radiation was filtered out 



by a series of dedicated algorithms 
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16| . They make use of differences between photon 



and cosmic-ray induced clusters in the topology-based signatures: the multiplicity of active 
modules in a cluster, the linearity of its shape, and the intrinsic dispersion of the energy dis- 
tribution in modules calculated with respect to the most energetic one. The invariant mass 
distribution of 77 pairs obtained is shown in the upper panel of Fig. [Hand exhibits a promi- 
nent 7r° peak. The invariant mass resolution (FWHM) of this peak is 11%, in agreement 
with Monte-Carlo simulations uM- The observed shoulder for masses below m^o is due to 



a partial escape of the electromagnetic cascade from a TAPS block 17(|. The high signal to 
background ratio S/B=18 allows selection of the 7r° candidates on an event-by-event basis. 

A kinematical fit [18] was applied to each event to improve the quality of the reconstruc- 
tion of neutral pion momentum from the momenta of the photons in a pair associated with 
the meson decay, taking advantage of the known neutral pion mass. The parameters of 
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this method, the invariant mass threshold m™" and the maximum chi-square of the event 
reconstruction Xmax^ were adjusted such that the distributions of pulls of the three mo- 
mentum components were in the best agreement with the Normal distribution. The best 
fit parameters were found to be m™" = 105 MeV/c^ and Xmax— 20. The influence of their 
variation on the results presented in this paper was added to the corresponding systematic 
errors. The bremsstrahlung photon pairs from the Ta+Au collision zone were automatically 
rejected by these conditions (mainly by the m™" cut). A total of 528 photon pairs satisfied 
all the above-mentioned filters and were assigned to vr^ decays. The contribution of cosmic 
radiation was found to be on the level of below one pair. 

The phase space population of measured neutral pions is shown in the bottom panel of 
Fig. [T]in the plane of rapidity and transverse momentum, where the midrapidity in the AA 
frame {Ucm) is indicated by an arrow. The dashed curve shows the average momentum of 
7r° emitted from a source obeying the Boltzmann distribution (Eq. [1]) with a temperature of 
15 MeV. This value is typical for the inverse slopes of momentum distributions of neutral 

nnn 

pions at beam energies around 40A MeV [Jj, QJ, |20l . The expected phase space distribution 



of neutral pions is well covered by the experimental data. To account for the detection 
efficiency, neutral pions were sampled from a Boltzmann distribution without collective 
radial flow of a particle |2l| 

-|^~p,.£.exp(-£/T), (1) 

where T is the temperature of the source, E = rriTC^ • cosh{y) is used and is the 



transverse mass mj-c^ = a^/ (ptc)'^ + {mrc^ |22i] . The response of the TAPS apparatus 



was simulated using the GEANT3-based j23l] KANE code 24]. The source temperatures 



were varied between 6 and 24 MeV, i.e. a range of inverse slopes of transverse momentum 
distributions found to be characteristic for beam energies between 25A and 95A MeV 
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27l |. The efficiency e of neutral pions was investigated as a function of their rapidity 
and transverse momentum (see the upper panel of Fig. [2]). In the region of px ^ 60 MeV/c 
the shape of e changes rapidly. Oscillations in rapidity stem from the geometric configuration 
of the TAPS blocks. In order to verify that the sharp drop of efficiency at lower transverse 
momenta is also due to the geometric setup of TAPS, a thermal source of 7r° mesons was 
generated within the PLUTO [28] code, and boosted according to the NN frame in the 
colliding system analysed. An efficiency Sgeom was constructed by selecting events, where 
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FIG. 2: Top: efficiency of vr'^ mesons as a function of rapidity and transverse momentum, obtained 
witli the KANE code [24i] . Bottom: estimation of geometric component of efficiency, simulated in 



frame of tlie PLUTO code 



281 ]. See text for details. 



both photons from a 7r° decay were emitted towards the TAPS blocks. The profile of 
as a function of pt, shown in the lower panel of Fig. |2l demonstrates that the TAPS 



-geom 



geometry is the driving cause of the drop observed in the upper panel of this figure. It 
also provides a cross-check of the GEANT-based analysis. In the regions of phase space 
around the limits of observed data some dependency of the efficiency on T was found. 
The existence of this dependency in conjunction with the profile of measured pt spectrum 
(c.f. Fig. |3]) is the reason behind using as an input to the KANE efficiency calculation a 
thermal distribution with a broad range of values of T parameter, rather than a homogenous 
spectrum m px — y phase space. Variations of parameters extracted in the present analysis 
due to the uncertainty of the efficiency map were added to the relevant systematic errors. As 
the efficiency changes quite sharply in some regions of pt and y, in order to obtain the weight 
for every reconstructed neutral pion, e was linearly interpolated between neighbouring bins. 

The overall 7r° acceptance of the apparatus and the applied analysis methods was found 
to be 0.36%. 
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IV. RESULTS 



In the previous measurement of Ta+Au collisions at 39. 5A MeV, performed by the TAPS 
collaboration, about 100 7r° mesons were reconstructed l|, l2|. The neutral pion cross section 
was found to be a^^o = 2.2±0.3 /ib and the transverse momentum spectrum was described 
by a Boltzmann function characterized by a temperature of 16±4 MeV. The cross section 
for bremsstrahlung photons of energies above 30 MeV, = 6.9±0.7 mb, was also obtained, 
thus the ratio of a^o/cr^ was found to be (3.2±0.5)-10~^. 

The analysis presented here is based on a number of identified pions 5 times larger 
compared to the previous measurement. Therefore, we are now able to study the properties 
of produced pions in more detail. As a cross-check with the previous measurement, the 
cTttO / cr^ ratio was obtained according to the following formula: 

^ ^ E.1'i i/g' 

where e'^^^ are the detection efficiencies of neutral pions and photons, and the summation 
runs over N^^/^ accepted particles. The relevant efficiency maps were simulated using the 
KANE code. The ratio obtained is found to be (4.5±0.4)-10~^, in 2a agreement with the 
value found in the previous measurement. For an absolute normalization of the neutral 
pion distributions presented here, the above mentioned cross section of 2.2±0.3 fib from 
Refs. [ll, l2| was used. 

For a narrow rapidity range Ay around Ucm, the Boltzmann distribution expressed by 
Eq. [1] can be approximated by: 

da ( rriTC^ {coshiy)) .\ 

^-PT-mr- exp |^ j , (3) 

To fit this function to the experimental data pions from the Ay = ±0.45 range were selected. 
This choice of rapidity range is a balance between the requirement of reasonable statistics 
and the reliability of approximation assumed in Eq. |3l Variations of the results presented 
below due to a choice of Ay were added to the systematic errors. The fit gives an inverse 
slope of 17.9±1.3j!'j''3 MeV (see solid curve in Fig. [3]^a) and Table [Tll)- However, a x'^/^ 
is found to be 2.5, and the auto-correlations in the residuals of neighbouring bins can be 
observed. The same discrepancy is also observed for the kinetic energy distibution, shown 
in Fig. |5]and discussed below. The fit of Eq. [3] to the high px part {pt> Q7MeV/c) delivers 





o 
> 

0) 

1. 
n 



Q. 



TJ 



, T = 17.9 ± 1.3 MeV 

T= 9.9 ±1.3 MeV 


- B 

/ 

. ,,,,,,, \ 


T = 21.8 

T = 12.3 


1 2.0 MeV 

1 2.1 MeV 



20 40 60 80 100 120 
Pt [MeV/c] 



FIG. 3: (Color online) Transverse momentum distribution of neutral pions emitted from Ta+Au 
collisions at a beam energy of 39. 5A MeV, for a Ay = ±0.45 window around , without (a) 
and including (b) the model of absorption in nuclear matter (see text). Thermal fits to the whole 
spectrum (solid lines) and to the high pT range (dotted lines) are shown. Arrows indicate the lower 
bound of the fit to a high px tail of the spectrum. 

an inverse slope parameter of T = 9.9±1.4l;Qg MeV, a value considerably lower than the 
one resulting from the fit in a full Pt range. However, in this case the low momentum 
part of the experimental spectrum is strongly overpredicted by the extrapolation of the fit. 
This depletion correlates with the decline of efficiency with p^-, as shown in the upper panel 
of Fig. [3l However, as reported in section IHIl an additional simulation verified that the 
positioning of TAPS blocks was solely responsible for the mentioned drop of efficiency. The 
findings reported above strongly suggest that the distribution of transverse momentum of 
neutral pions cannot be fitted globally with the Boltzmann's thermal function. As shown 
in Table [H this trend has been observed in a ran ge of sy stems investigated by the TAPS 
collaboration at beam energies below 100 A MeV jll. li9l.l2O. 25-27|. The pt spectra of neutral 
pions from other experiments were taken as they were published and fitted both to the full 
range and with a common lower bound of 60 MeV/c with the function 



/(Pt) ~ Pt ■ ■ exp{-mT/T) 



(4) 



In spite of somewhat different ranges of rapidities accessible to those experiments, in all 
reported cases the slope obtained globally was larger than the one fitted to the high pt tail 
of distribution. Therefore, as for the current data, extrapolations of Boltzmann functions 
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TABLE I: Systematics of inverse slopes obtained by fitting tlie Eq. |3]to the transverse momentum 
distributions of neutral pions at beam energies below lOOA MeV from recent experiments. Global 
fits {pt > MeV/c) and high pT (> 60 MeV/c) are shown. In the cases where authors used 
different functions, or lower fit bounds for fitting the tail was different than 60 MeV/c, the data 
were extracted and refitted. 



System 


Tbeam [A MeV] 


T (pT > MeV/c) [MeV] 


T {pT > 60 MeV/c) [MeV] 


Reference 


Ar+Au 


25 


13.3+2.5 


6.4+1.9 


[26], prelim. 


Ar+Au 


35 


13.0+0.9 


11.6+0.9 


[26], prelim. 


Ta+Au 


39.5 


17.4+1.31?;^ 


11.0+1.2li;0 


this work 


Xe+Au 


44 


14.6+1.0 


12.7+1.0 


[19] 


Kr+Ni 


60 


17.0+0.5 


15.8+0.2 


[1] 


Ar+C 


60 


16.2+0.7 


13.3+1.0 


[25, 3^ 


Ar+Ni 


60 


16.6+0.6 


11.9+0.3 


[25, 3^ 


Ar+Ag 


60 


18.3+2.8 


16.3+2.7 


[25,31] 


Ar+Au 


60 


17.2+0.8 


12.7+0.7 


[25, 31] 


Ar+C 


95 


16.1+0.2 


14.2+0.1 


[25, 3^ 


Ar+Al 


95 


18.4+0.1 


16.4+0.1 


[25, 3JJ 


Ar+Ag 


95 


16.8+0.2 


14.1+0.2 


[25, 3JJ 


Ar+Au 


95 


15.2+0.2 


13.6+0.2 


[25, 3JJ 



fitted to the high pt tails overpredicted the measured pion yield at lower transverse momenta. 
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30|, no depletion of 



In the measurements performed in the eighties at 25A and 35A MeV 
yield was found. However, these experiments were performed with lead-glass spectrometers, 
which allowed for an invariant mass resolution of only 40%. Also, the experimental setup 
used in these measurements did not allow for the application of advanced procedures for the 
rejection of cosmic-ray background, employed in the analysis presented in our paper. 

To estimate the effect of the shadowing of pions by nuclear matter, a simple geometrical 
model 
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33| was applied. It assumes 7r° production from a random point inside a static 



collision zone defined by the maximum overlap of two spheres at a random impact parameter. 
Mesons are produced in a thermalized source of temperature Tmodei = 12 MeV and emitted 
with an angular distribution: 
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FIG. 4: Experimental angular distribution of neutral pions emitted from Ta+Au at 39. 5A MeV 
shown in the NN frame, compared to Xe+Au at 44A MeV 
the second order Legendre polynomial (see text). 



19l |. Curves represent the best fits of 



f{^) = l + A2-P2{cos^), (5) 

where P2 is the second order Legendre polynomial and A2 = 0.05 is the anisotropy 
parameter 35| . This model uses the momentum-dependent mean free path of a neutral 
pion js^. However, it must be noted that neither pion scattering nor the dynamics of the 
collision is included in this model. 

The primordial pt distribution estimated by this model is shown in Fig. [3](b) . A fit of 
Eq. |3] performed in the whole range of transverse momenta delivers the inverse slope of 
21.8±2.0 MeV and a x^/^ = 1-6- This is an improvement compared to x^/'^ = 2.5 for the 
measured data, so an application of this model brings the distribution somewhat closer to the 
thermal shape. However, auto-correlations in the residuals of neighbouring bins still persist. 
Fitting the high pt tail with a Boltzmann distribution gives a considerably lower inverse 
slope of T = 12.3±2.ltJ i MeV. However, the fitted function still strongly overestimates 
the measured yield at lower pt values. In the scenario of thermal equilibrium established 
during the collision, and the system cooling gradually, neutral pions of high pt or kinetic 
energies should be rather characterized by highest temperatures. However, this approach 
is in contradiction with the findings presented above. The mismatches demonstrated above 
suggest that neutral pions produced and emitted from the collision zone may not be ther- 
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mally equilibrated, even if the absorption of neutral pions is accounted for in frame of the 
above mentioned model. 

The experimental polar angle distribution of neutral pions in the NN frame, (see Fig. H]) 
exhibits an anisotropy with less preference for emission in sideward directions. The lack 
of strong asymmetry between forward and backward directions can be explained by the 
near-symmetry of the colliding system. The pronounced sideward absorption can be quali- 
tatively explained by the longer path of a pion in the colliding nuclear matter due to mutual 
shadowing by the heavy collision partners. The fit of Eq. |5] to the experimental data yields 
an anisotropy parameter A2 = 0.31±0.09!!^og2- In an attempt to reconstruct the primordial 
polar angle distribution of neutral pions, the above-mentioned 7r° absorption model was 
used. To parametrize the primordial spectrum, a dipolar distribution of the form of Eq. \5\ 
was assumed, with an anisotropy parameter A^^^. The angular spectrum of pions which 
traversed the hadronic matter was compared to the spectrum observed by TAPS. The best 
fit to the experimental data was obtained for A^^"^ = 0.09±0.09^q;q2, which exhibits a near 
isotropic emission pattern. Despite the meager experimental knowledge of ^^""^ values for 
collisions of symmetric systems at comparable beam energies, an average anisotropy param- 
eter AiF^^'^ = 0.32±0.05 for the primordial distribution of all hitherto investigated systems 



was found 
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33| . This value is quite higher than the one obtained for Ta+Au. However, 



several values are also beyond the range of the global average quoted above. A fit of Eq. [5] to 
ihe angular distribution of neutral pions emitted from ^^^Xe+^^^Au collisions at 44A MeV 
19| (see bottom panel of Fig. H]) provides an anisotropy parameter of A2 = 0.30±0.04. An 
estimation of this parameter for the primordial pions done in the framework of the above 
mentioned 7r° absorption model gives Ag"™ = 0.10±0.04. Both values are in very good 
agreement with results for Ta+Au collisions obtained in this work. 

The kinetic energy distribution in the NN frame has also been investigated. Global 
fitting by a Maxwellian distribution of the form ~ pE ■ exp{—E/T) [see Fig. [5t^a)] results 
qualitatively in the same picture as in the case of px spectrum: high value of x^/^ = 4.5 
and auto-correlations in the residuals of neighbouring bins. A fit of the same function to the 
high energy tail (i?fcj„>37 MeV) gives an inverse slope of T = 8.7±1.6;!^53 MeV. Similarly 
to the case of the pt distribution, the application of the above mentioned vr'' absorption 
model results in a slightly increased inverse slope of the high energy tail of the estimated 
primordial spectrum, T = 10. 8±2. 215^7 MeV, and the depletion of the yield in the low E^in 
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TABLE II: Summary of inverse slopes and polar anisotropy parameters. Also, shown are the estimates of the primordial values of these 
parameters, calculated within the 7r° absorption model described in the text. For the contributions to the systematic errors, see text. The 
global fits {pt > 0, Ekm > 0) could not describe the spectra well. 
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FIG. 5: (Color online) Kinetic energy spectrum of neutral pions, without (upper panel) and in- 
cluding (lower panel) the correction based on the vr" absorption model (see text). Curves represent 
fits of Maxwellian distribution to the whole spectrum (solid lines) and to the high -Efcm. part (dot- 
ted lines; vertical lines mark the lower fit bounds. Arrows indicate the kinematical limit for the 
NN —7- NNtt^ reaction including Fermi motion. The TAPS response function for pions emitted 
with this limiting energy is shown by the dashed lines. 

part of the spectrum cannot be removed [see Fig. Et^b)]. 

The constructive superposition of the beam momentum and two Fermi momenta of nu- 
cleons {pf = 270 MeV/c) from the colliding nuclei enables the production of pions up to 
kinetic energy of E^^^^ = 37 MeV. However, the measured kinetic energy spectrum of neutral 
pions is spread far beyond this value. The experimental smearing of this upper energy limit, 
simulated using the KANE code 2J] , and shown in Fig. [5l cannot explain the extent of the 
high Ekin shoulder. One possible hypothesis leading to an explanation of this effect could be 
the production of pions in channels different from the process A^A^ — )■ NNtt^, like the multi- 
step process i) A^A^ A^A, AA^ NNir^ or the cooperative one ii) A^A^A^ NNNtt'^. To 
estimate the kinetic energy limit for each of these processes, which would also take into ac- 
count the momentum conservation of the reaction ^oducts involved, the collisions studied 
were simulated with the help of the PLUTO code [28]. A kinetic energy of 39. 5A MeV was 
assigned to a "beam nucleon". The Fermi motion of nucleons was included. In the second 
step of channel i) all the available energy was assumed to be transferred to the A(1232) 
resonance mass, an assumption justified by the rise of the A production probability with a 
resonance mass in the region well below {mj\). The A resonance interacted with another 
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nucleoli of Fermi gas, originating randomly from the target or beam nuclei. We found that 
the kinetic energy limits amount to around 75 MeV for both investigated channels. This 
may be a hint that the investigated channels could play a role in explaining the observed 
presence of neutral mesons above the -E^^ limit. On the other hand, calculations within 
framework of the Dubna Cascade Model, assuming vr^ production in NN collisions, and 
employing inelastic vrN scattering {sj], systematically underestimated the measured yield in 
the high energy tail of the spectrum. More systematic investigations of possible dynamical 
processes are needed to explain the experimental findings. 

V. CONCLUSIONS 

Neutral pions produced in the Ta+Au collision zone at 39. 5A MeV were measured using 
the TAPS electromagnetic calorimeter. Their transverse momentum distribution selected 
around midrapidity and kinetic energy spectrum deviate from a thermal shape. These devia- 
tions cannot be removed by the application of a phenomenological model of 7r° absorption in 
nuclear matter. The kinetic energy spectrum reaches far higher values than the limit for the 
NN — )■ NNtc^ channel, including constructive boost of Fermi motion. It is suggested that 
multi-step or cooperative processes could be possible sources of the neutral pion production. 
Nonetheless, the tails of the pt and E^in spectra were in agreement with Boltzmann distribu- 
tions, yielding inverse slope parameters of T = 9.9±1.4+J:° MeV and T = 8.7±1.6+°:^ MeV, 
respectively. The polar angle distribution of the measured pions reveals a dipolar anisotropy 
characterized by A2=0.31±0.091o'q2- However, the estimated asymmetry of the primordial 
angular distribution, calculated within the 7r° absorption model, turned out to be almost 
isotropic (^r™ = 0.09±0.09l°:°^). 
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